Weak momentum scattering and the conductivity of graphene 
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Electrical transport in graphene offers a fascinating parallel to spin transport in semiconductors 
including the spin-Hall effect. In the weak momentum scattering regime the steady-state density 
matrix contains two contributions, one linear in the carrier number density n and characteristic 
scattering time r, the other independent of either. In this paper we take the Liouville equation as 
our starting point and demonstrate that these two contributions can be identified with pseudospin 
conservation and non-conservation respectively, and are connected in a non-trivial manner by scat- 
tering processes. The scattering term has a distinct form, which is peculiar to graphene and has 
important consequences in transport. The contribution linear in r is analogous to the part of the 
spin density matrix which yields a steady state spin density, while the contribution independent 
of r, is analogous to the part of the spin density matrix which yields a steady state spin current. 
Unlike in systems with spin-orbit interactions, the n and r-independent part of the conductivity 
is reinforced in the weak momentum scattering regime by scattering between the conserved and 
non-conserved pseudospin distributions. 
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I. INTRODUCTION 

The zero-gap semiconductor graphene, or two- 
dimensional carbon, is a new material with a host of re- 
markable physical properties that offers the possibility of 
all-carbon devices. The last three years have seen a surge 
of experimental and theoretical interest in graphene fol- 
lowing its realization in the laboratory. ^"■^ High- mobility 
graphene samples are nowadays reliably manufactured, 
and the recent experimental success in applying a top 
gate* offers an increased handle on material properties. 
Among the latest notable successes the ability to suspend 
graphene experimentally"" is expected to help determine 
the value of the dielectric constant e unambiguously in 
this material. 

The linear spectrum characterizing the band structure 
of single-layer graphene is reminiscent of dispersion rela- 
tions in relativistic physics, and the charge carriers in this 
material behave like massless Dirac particles. Graphene 
has a honeycomb lattice with two atoms per unit cell, 
and the Hamiltonian displays a coupling between elec- 
tron and hole states which gives rise to a degree of free- 
dom we refer to as the pseudospin. It is these facts that 
underlie its unusual features, which include the vanishing 
of the density of states at the Dirac point, a contribution 
to the conductivity independent of the carrier density n 
and scattering time r, a half-integer quantum Hall effect, 
and the Klein paradox. ^"'^ 

Recent experimental work includes the fabrication 
of epitaxial graphene/graphene-oxide junction,® mea- 
surement of ultrafast carrier dynamics,^ shot noise 
measurements,*'® determination of the performance lim- 
its of graphene devices, ""^^ observation of the Aharonov- 
Bohm effect, ^^ observation of the quantum Hall effect 
near the Dirac point, -'^^ of a rcnormalization of the ve- 
locity due to electron-phonon interaction,^^ and a thor- 
ough experimental study of epitaxial graphene on SiC.^* 
Theoretical research on single-layer graphene has concen- 



trated on, among other matters, the effect of electron- 
electron interactions,^'^ which are expected to be impor- 
tant due to weak screening, on the question of whether 
graphene is a Fermi liquid,^® and on the importance of lo- 
calization around impurities. ^^ A number of theories have 
dealt with scaling, ^^'^^ impurity states, ^° the odd-integer 
quantum Hall effect, ^^ the fractional quantum Hall 
effect, •^^'^•^ polaritons,^* spin-orbit coupling^^ and sum 
rules for the optical and Hall conductivities.^® In addition 
to these, a large number of theoretical predictions include 
the quantum spin-Hall effect, ^^ spin-Hall conductance 
fluctuations,^* proximity-induced superconductivity,^® 
antiferromagnetism,^" the spin-valve effect, ^^ peculiar 
focusing properties of graphene p-n junctions, '^^ the 
use of graphene quantum dots as spin qubits,'^'^ Weiss 
oscillations,'^ and a zero-bias anomaly in the tunneling 
density of states. ^^ Many other theories have sought to 
increase theoretical understanding of graphene. ■^®"®'^ Be- 
yond the single-layer form, graphene ribbons have been 
predicted to have spin polarized edge states,®* while the 
band gap in bilayers has been shown to be tunable by 
means of an electric field. ®^ 

In the following we shall refer frequently to the weak 
momentum scattering regime, characterized by sp 3> 
fi/r, where r is a characteristic momentum scattering 
time, and the strong momentum scattering regime in 
which ep ^ h/r. Furthermore, it is conventional in the 
literature to make the distinction between intrinsic and 
extrinsic graphene. Intrinsic graphene refers to the spe- 
cific case in which the carrier doping density n is zero, 
and the Fermi energy lies at the Dirac point, k = Q. In- 
trinsic graphene is therefore by definition in the strong 
momentum scattering regime. Extrinsic graphene refers 
to the doped case and may be in either the weak or the 
strong momentum scattering regime. 

The presence of the additional degree of freedom con- 
tained in the pseudospin causes the steady state for 
graphene in an electric field to be qualitatively different 



from any other known material. In partieular, the pres- 
ence of a contribution to the conductivity independent of 
the carrier density remains a puzzhng observation. Such 
a contribution has been measured experimentaUy^"'^ and 
we emphasize that it was extracted from the experimen- 
tal data by taking the doping density to zero. The num- 
ber obtained is naturally characteristic of the strong mo- 
mentum scattering regime. It is typically referred to as 
the minimum conductivity of graphene. At the same 
time, theoretical research on clean graphene (no scatter- 
ing) has found an additional contribution to the conduc- 
tivity independent of n and t.^^^^^ There appears to be 
some agreement that coherence between electrons and 
holes lies at the heart of this particular contribution to 
the conductivity. Since the Hamiltonian for carriers in 
graphene is hvcr ■ k, where w is a constant and cr repre- 
sents the pseudospin, the velocity operator is simply vcr, 
and it is evident that the pseudospin will play a crucial 
role in transport. Nevertheless this contribution to the 
conductivity was found in the ballistic regime, whereas 
the contribution referred to as the minimum conductivity 
was measured experimentally in the strong momentum 
scattering regime. We note that a series of enlighten- 
ing papers have focused on Boltzmann transport^^, on 
transport in the strong momentum scattering limit^^"'*^ 
and transport in extrinsic graphene*^'^^. Our work differs 
from these papers in that it is the first to recover both the 
ordinary r-dependent "Boltzmann" conductivity and the 
n- and r-independent conductivity, and to demonstrate 
the profound relationship that exists between the two. In 
addition, our formalism does not require us to work in the 
limit n — > in order to recover the n- and r-independent 
conductivity. 

In order to arrive at a tractable equation for the density 
matrix it is necessary to assume that the Fermi energy 
ep 3> ft/r, which requires us to restrict our attention 
to the weak momentum scattering regime. We assume 
low temperatures, where scattering due to charged im- 
purities is important and may be dominant, and the fact 
that weakly doped graphene is a Fermi liquid^^ justifies 
omitting the effect of electron-electron scattering. We 
believe that our approach sheds light on the manner in 
which the part of the conductivity of graphene arises that 
is independent of number density and t and how it is re- 
lated to the ordinary, r-dependent ("Boltzmann") con- 
ductivity. Indeed, the theory presented in this work is 
tailored towards rendering explicit the important role of 
the pseudospin in the dynamics of carriers in graphene, 
with a focus on steady state processes, providing an ac- 
curate and at the same time transparent approach. We 
wish to emphasize that, although we determine the nu- 
merical value of the conductivity, our principal aim is not 
to obtain a number but rather to bring to light the un- 
derlying structure of the steady-state density matrix in 
pseudospin space. Compared to the Boltzmann picture, 
our formalism has coherence between electrons and holes 
built in from the start and thus provides a clear phys- 
ical picture. Nevertheless weak localization effects'*®''*^ 



are not taken into account in this work, as they are not 
expected to be important in the weak momentum scat- 
tering regime. 

We find the most important observation to be the fact 
that a carrier's pseudospin is not conserved, because of 
electron-hole coherence present in the Hamiltonian (in- 
cluding coherence induced by the electric field). In the 
absence of intervalley scattering, assumed in our work, 
a pseudospin eigenstate is an electron or a hole, thus 
pseudospin non-conservation means a continually chang- 
ing combination of an electron and a hole. Therefore each 
carrier can be thought of as a part which is either an elec- 
tron or a hole, and a part which is a continually changing 
mixture of an electron and a hole. With this in mind, it 
makes sense to divide the pseudospin density matrix into 
two linearly independent contributions, corresponding to 
conserved pseudospin (the carrier is an electron or a hole) 
and non-conserved precessing pseudospin (the carrier is a 
continually changing mixture of an electron and a hole). 
We show in a systematic fashion that the two contri- 
butions to the conductivity of single-layer graphene are 
related to these linearly independent components of the 
pseudospin density matrix. The r-dependent contribu- 
tion is a result of pseudospin conservation, while the n- 
and r-independent conductivity stems from pseudospin 
non-conservation. The two independent parts of the den- 
sity matrix, often referred to as the dissipative and re- 
active parts, which are responsible for the two terms in 
the conductivity, are connected in a non-intuitive way 
by scattering events. We find that scattering gives rise 
to a term in the Liouville equation which is peculiar to 
graphene and distinct from the usual scattering term in 
other conductors. In particular, even an infinitesimal 
amount of scattering produces a contribution to zeroth 
order in the scattering potential which reinforces the n- 
and r-independent contribution to the conductivity com- 
ing from the band structure. This feature is one of the 
many intriguing analogies which exist between transport 
in single-layer graphene and the generation of steady- 
state spin densities and currents, including the spin-Hall 
effect, in ordinary semiconductors. ^^'''^'^^ From a differ- 
ent perspective, the two independent parts of the density 
matrix are intimately related with the phenomenon of 
zitterbewegung which refers to a highly oscillatory com- 
ponent in the motion of rclativistic particles. ^^ Quite 
generally, such an oscillatory component is always ob- 
served when the electron states in two or more neigh- 
boring bands interfere.''^ Of course, such an interference 
also hes at the heart of (psuedo)spin precession so that 
the unique transport properties of graphene can be con- 
sidered to be a manifestation of zitterbewegung in a solid- 
state system.**^ 

The outline of the paper is as follows. In Sec. II we 
derive a kinetic equation for the density matrix specific 
to graphene, taking the quantum Liouville equation as 
our starting point and focusing on the form of the scat- 
tering term, which is different in graphene from that in 
other materials even in the first Born approximation. Fol- 



lowing that, in Sec. Ill we consider the dynamics of the 
pseudospin. We divide the pseudospin density matrix 
into a part representing conserved pseudospin and a part 
representing non-conserved pseudospin, and demonstrate 
the way scattering affects both of these parts and con- 
nects one to the other. In Sec. IV we determine the 
steady-state solution in the presence of an electric field 
and the electrical conductivity. We show that the ordi- 
nary, T-dependent conductivity can be traced to the con- 
served pseudospin distribution, while the number den- 
sity and T-independent conductivity is associated with 
the non-conserved pseudospin distribution. These argu- 
ments are further developed in Sec. V. We illustrate in 
Sec. VI the remarkable similarities between charge trans- 
port in graphene and spin transport in semiconductors 
with strong spin-orbit interactions, including the way 
the vertex correction to spin currents in semiconductors 
has an analog in graphene that reinforces the n- and r- 
independent conductivity in the weak momentum scat- 
tering regime. We conclude with a brief summary in 
Sec. VII. 



II. TIME EVOLUTION OF THE DENSITY 
MATRIX 

The system is described by a density operator p which 
obeys the quantum Liouville equation 



f + ^[^ + ^" + ^'^]-0' 



(1) 



where H is the band Hamiltonian, H^ represents the 
interaction with external fields, and U is the impurity 
potential. We project the Liouville equation onto a 
set of time-independent states of definite wave vector 
{|fcs)} that are not assumed to be eigenstates of H. 
The matrix elements of p in this basis are written as 
Pkk' = Pkk' — {^s\pWs')j with corresponding nota- 
tions for the matrix elements of iJ, H^ and t/, thus 



Pkk' 



Hi 



kk' 



Huu,, and Ukk' are 2x2 matrices in the 



space spanned by the pseudospin s. We refer to pkk' as 
the (pseudospin) density matrix. Matrix elements of the 

^kk' 



band Hamiltonian Hkk> = Hk S^k' and H^u, = HP S^k' 



are diagonal in k but contain off-diagonal terms in the 
pseudospin indices. Matrix elements of the scattering 
potential Ukk' are off-diagonal in k. (Matrix elements 
of the form Ukk lead to a redefinition of Hk-) We as- 
sume elastic scattering and work in the first Born ap- 
proximation, in which C/^^/ = Ukk'Sss'- Impurities are 
assumed uncorrelated and the normalization is such that 
the configurational average of {ks\U\k's'){k's'\U\ks) is 
{ni\Ukk'\'^Sss')/V , where n^ is the impurity density, V 
the crystal volume and U^k' the matrix element of the 
potential of a single impurity. 

Pkk' is divided into a part diagonal in k and a part off- 
diagonal in fe, given by pkk' = fkSkk' + 9kk' , where, in 
9kk' : it is understood that k ^ k' . We will be interested 
primarily in fk since most operators related with steady 
state processes are diagonal in k. The quantum Liouville 
equation is broken down into 



'^ + \[HkJk] = -Hi-^[U,g]kk, (2a) 



dg, 



kk' 



^^liH,g]kk' = -^iUj + gU' 



(2b) 



The solution to Eq. (2b) to first order in U can be written 
as 



^ Jo 



1,1 —iflt' Ih 

gkk' = - I at e ' 



uJit-t') 



e ' Ifcfc' 



(3) 
The assumption that Sft/H ^ 1 allows us to expand 

f{t — t') around t and retain only /(t). (Additional terms 

are of higher order in U.) The equation for //, is 



^ + l[Hk,fk] + J{fk) = -Hi, 



1 r°° 
Jifk) = j^l dt' 



U,e 



-iHt'/h 



uj{t) 



iHt'/h 



(4a) 



i*' 



The integral in Eq. (4b) is performed by inserting a reg- 
ularizing factor e~^* and letting 77 — > in the end. For 
potentials | Ukk' \ oc 1 which are scalars in pseudospin 
space this integral has the form 



J 



Jifk) 






^ \Ukk'\' I dt' e'^^' [e-^"-'' /^{,fk - f'k) e'""'^^ + e'^"-' '^{fk - fL) e^"-''/^}. 



(fk' 



(5) 



Equation (4) is a generalization of Fermi's golden rule or, 
equivalently, a generalization of the first Born approxi- 
mation to systems, where the orbital motion is coupled 
with a (pseudo)spin degree of freedom. ®^"^" 



III. PSEUDOSPIN DYNAMICS 



In the following we derive a scattering term specific to 
graphene and an equation describing the time evolution 
of the pseudospin. The band Hamiltonian for the carri- 



ers in each valley in single-layer graphene at low doping 
densities is given by H^ = hv cr ■ k, where the constant v 
is the Fermi velocity and cr is the (two-dimensional) vec- 
tor of Pauli matrices in pseudospin space. We emphasize 
that the Hamiltonian docs not depend on the true spin 
of particles, thus the final result will contain a factor of 
2 from the sum over the spin. An additional factor of 
2 must account for the valley degeneracy. Consequently 
final results are multiplied by an overall factor of 4. 

The Hamiltonian H^ is formally similar to the spin- 
orbit interaction in spin-1/2 electron systems, ^^ except 
that the spin-orbit interaction is usually accompanied by 
a kinetic energy term quadratic in k, which is typically 
much larger than it, and has no analog in graphene. We 
wish to consider briefly this aspect of Hk from the point 
of view of symmetry. The Hamiltonian Hp, transforms 
as a dipole in pseudospin space,® unlike the Hamiltoni- 
ans of spin-1/2 electron systems, which contain both a 
dipole (the spin-orbit interaction) and a monopolc (the 
scalar kinetic energy). In the absence of scattering the 
equations of motion for the monopole and the dipole are 
decoupled,®^ but when scattering is present the time evo- 
lution operator in the scattering term (5) in general mixes 
the monopole and the dipole. In graphene the monopole 
of the density matrix is equivalent to the scalar part and 
the dipole is equivalent to the pseudospin part. Inter- 
estingly, because the Hamiltonian for graphene only con- 
tains the dipole term, scattering cannot mix the scalar 
and pseudospin parts of the density matrix. We will see 
in the next paragraph how these symmetry arguments 
become relevant. 



A. Scattering term 

Next we will evaluate the scattering term (5). For this 
purpose, the density matrix of graphene is written as 
fk = nk + Sk = rikt + ^ Sk ■ c with a scalar part rik 
and pseudospin part Sk- We substitute H^ in the time 
evolution operator in Eq. (5) and carry out the time inte- 
gration. In this way we obtain expressions for the action 
of J on rifc and Sk, and we flnd that, as discussed above, 
this term does not mix nk and Sk- The explicit expres- 
sion for J{nk) is not needed in this work and will not be 
given. The action of the scattering term on Sk is 



J{Sk) 



I de'\Ukk'?{Sk-Sk' 

Jk'=k 



krii 
Shirv Jk'=k 

[a {1 - cosj) + {(T - k) k' + k{(T-k')]- (6) 



In the above 0' is the polar angle for the direction of k' 
and 7 = 6*' — 6* is the angle between k' and fc. We see in 
Eq. (6) that the scattering term is qualitatively different 
from the scattering term in spin-1/2 electron systems,®^ 
both in its angular dependence and in not mixing the 
scalar and pseudospin parts of the density matrix. Aside 
from the Born approximation, no further approximations 
were made in deriving Eq. (6). Yet we emphasize that 



it is essential for this theory to assume weak momentum 
scattering since there is no equivalent of the scalar kinetic 
energy term present in semiconductors with spin-orbit 
interactions, and in order to derive Eq. (4) one must as- 
sume EpT /h = vkpT ^ 1. The problem is characterized 
by two time scales, the pseudospin precession frequency 
ep/fi = vkp and r, and in assuming vkpT 3> 1 we were 
able to truncate the scattering term at the leading order 
in the impurity potential. 



B. Time evolution of the pseudospin 

We consider in more detail the equation for the pseu- 
dospin part of the density matrix Sk, in the general case 
in which a nonzero source term exists on the RHS. Such a 
source term will be present when an external field is act- 
ing on the system, and its form can be derived straight- 
forwardly from the quantum Liouville equation. The spe- 
cific case of an electric field will be discussed in the follow- 
ing section. The kinetic equation for Sk in the presence 
of a source S^ is 



dSk 
"dT 



[H-k, Sk\ + J{Sk) 



(7) 



The structure of this equation is very important. In or- 
der to bring it out we decompose Sk into two linearly 
independent parts, Sk = Sk\\ + Sk^- Sk\\ is, in matrix 
language, parallel to the Hamiltonian Hk, while the re- 
mainder Ski_ is orthogonal to Hk- Because [Hk, Sk\\\ = 
the parallel part Sk\\ does not change in time under the 
action of the time evolution operator e^^kt/h^ jj^ other 
words Sk\\ represents the fraction of the pseudospin which 
is conserved, that is, carriers which are either electrons 
or holes. Conversely, Sk± represents the fraction of the 
pseudospin which is not conserved, i.e., it is precessing. 
This fraction corresponds to carriers which are a continu- 
ally changing mixture of electrons and holes. The decom- 
posotion reficcts the central importancee of electron- hole 
coherence for the carrier dynamics in graphene. We re- 
mark that the decomposition is fully equivalent to the de- 
composition of carrier dynamics in relativistic quantum 
mechanics which contains a smooth part and on oscilla- 
tory part known as zitterhewegung.'^'^'''^ In that sense the 
physics discussed here represents a direct manifestation 
of zitterbewegung in graphene.'*'^ 

Specifically, the following equations hold for Sk\\ and 
Sk± 



Skii — 



tr(5fci7fe) 



triH'k 



Hk 



[Sk\\,Hk] 
tiiSk^Hk) 



0, 
0, 



(8a) 

(8b) 
(8c) 



where the symbol tr refers to a trace over the pseudospin 



indices only. It is easily seen that 

Sk\\ = 5 (S'fc • fc)(cr -fe) EE isfe||crfc||, (9a) 

Sfcx - ^{Sk-e){(7-e) = ^sk^ak±, (9b) 

where k and 9 are unit vectors along the direction of k 
and perpendicular to k, respectively. We note that any 
matrix in pseudospin space can be decomposed as in Eq. 
(9). Therefore, in analogy with this decomposition of the 
pseudospin part of the density matrix, we also decompose 
the source term S^ along the same principles. From Eq. 
(7) we can immediately see that the equations describing 
the time evolution of S'^ii and Sk± are 



dt 



+ P\lJ{Sk) = Sfcii, (10a) 



dSi 



k± 



dt 



[Hk,Sk±]+P±JiSk) 



(10b) 



where Pn/^J(Sk) indicates that the scattering term acts 
on Sk = Sf,\\ + Ski_ and the resulting expression is pro- 
jected parallel/perpendicular to the Hamiltonian. 

A solution for Sk can be found most straightforwardly 
by expanding S^w and Sk± in the transition rate \Ukk'\'^, 
in a manner analogous to that adopted in determining the 
steady states of spin distributions in systems with spin- 
orbit coupling. We found^^ that in steady-state prob- 
lems the density matrix always contains a correction oc r 
and is thus of order —1 in the transition rate. This tells 
us that the expansion of Sk needs to start at order —1. 
Since we are working in the weak momentum scattering 
limit we truncate this expansion at the next highest or- 
der, which is order zero. The source term Sfc does not 
have any dependence on the transition rate and is thus 
of order zero. Equating terms of the same order in the 
transition rate in Eq. (10) shows that the expansion of 
Sf,\\ must start at order —1, while the expansion of Ski_ 
must start at order zero. We denote the order —1 in 
the transition rate by a superscript (—1) with the corre- 
sponding notation for order zero. As a result Eqs. (10) 
in the weak momentum scattering limit simplify to 



^11 ^(^i'^) 



^k\\, 



(11a) 



dS 



(0) 
k± 



dt 



[Hk,S'^ 



(0)l 



P|l J{S, 



(0)^ 
fell '' 



^k±-P±JiSk\\) (lib) 
-P|l J(5S). (lie) 



In Eq. (11a) we have omitted the time-derivative of S^n 
for the following reason. In both equations we are look- 
ing for the steady state solution, and the equation for 
Sf,\\ is most easily solved without the time derivative. 
The equation for Ski_ is most easily solved with the time 
derivative explicitly taken into account, but the time de- 
pendence drops out in the end. 

Equation (11) shows that, if the solution is required 
only to order zero in the transition rate, the scattering 



term acts only on Sk\\, the part of the density matrix 
parallel to the Hamiltonian. Physically, the fact that 
Sk\\ starts at a lower order in the transition rate than 
Sk± means that scattering processes are more effective 
at randomizing the pseudospin than at scattering into 
pseudospin eigenstates. 

Finally, we can simplify the scattering term by project- 
ing Eq. (6) onto and perpendicular to the Hamiltonian 
Hk- The projections that we will require in this work are 



P\\JiSk\\) 
P±JiSk\\) 
P\\J{Sk±) 



krii 
Shnv, 

kn. 
Shnv, 

kn. 



de' |C/fefe'P (sfcll - Sfc||)(l + cos7)crfc||, 

(12a) 

■ / d9' \Ukk'? (sfeii - 4||) sinTc^feX: 



(12b) 



— -^ d0' \Ukk'\'^ {sk± + Sk'±) sin7crfc||. 



(12c) 



We proceed to determine the concrete form oi J(Sk\\) for 
a screened Coulomb potential. The effect of screening in 
graphene has been evaluated by Ando''^ among others, 
who showed that k^p oc kp, where /ctf is the Thomas- 
Fermi wave vector. The explicit expression kxp will not 
be reproduced here, it suffices to bear in mind that the 
ration kxp/kp is a constant. In two dimensions, the 
square of the matrix element Ukk' of a screened Coulomb 
potential between plane waves is 



It/, 



1 



fcfc'l — 



elVi^Ak^.in^ 
W 



1/i? 



sin2(7/2) + fc2 



Tf/ '^F 



(13a) 
(13b) 



where Z = 1 is the ionic charge, Ls — kp / (2kTpk) is 
the screening length. Substituting this into Eq. (12a) we 
obtain 



P\\JiSk\\) 
C(7) 



kn^W 

cos2(7/2) 
sin2(7/2) + fc2p/fc| 



de'ai)isk\\-s'k«)cTk\\ (14a) 



(14b) 



with similar results for P^J{Sk\\). In order to evalu- 
ate this expression we expand C(7) in a- Fourier series as 
C(7) = SmCme™'^ and remark that C-m = Cm- In a 
similar way we expand s^n as s^n = ^^ Sfc||m e*"*^. This 
gives for Eq. (14) 



P\\J{Sk\\ 



kriiW 



7 ,(Co ^ Cm) Sk\\r 



Am9 



(15) 



This is the furthest this equation can be simplified at 
this stage. We will see below that in the steady state 
additional simplifications emerge. 



IV. STEADY STATE SOLUTION 



B. Sk± to leading order in scattering 



In the following, we assume low fields E and look for 
a solution to first order in E. As shown in Appendix A 
in the presence of an electric field E the source term S^ 
in Eq. (7) takes the form 



Sfc — — — • 



eE dSo 



dk 



(16) 



Here So is the pseudospin part of the equilibrium density 
matrix, i.e.. 



•So — 2 (/o+ ^ /o-) CTfell 5 



(17) 



where the scalars /o± = fo{±hvk), with /o the Fermi- 
Dirac distribution and ±hvk the eigenenergies of the 
graphene Hamiltonian Hf,, thus 



/o± — 



1 



gf3(^hvk-fi) j^ I ' 



(18) 



where /i is the chemical potential. The conserved and 
non-conserved components of the source term are 



eE-k fdfo+ dfn- 



2ft V dk 



dk 



ffcll 






A. Sk\\ to leading order in scattering 



(19a) 
(19b) 



^(0) 



We proceed to determine >S'^j_. For this purpose we 

fell 



require the term P± J(>S'L ) in Eq. (lib). Inserting Eq, 



(22) into Eq. (12b) gives 



(_i) ^ eE-e ^0 f dfo+ _ dfo- 
II ' 2ft Co - Ci V dk dk 



where the angular integral ^o is given by 

de' 



ffei , 

(23) 



6 



1 



2tt 



■ 2 -V 2 "v 

sm ^ cos -J 



n2 7 




fc|p fcxF 



(24) 

so that ^o/(Co ^ Ci) = 1/2 and Eq. (23) is indeed ind- 
pendent of the screening length. Now Eq. (lib) for Sf,j_ 
becomes 



dS: 



(0) 



eE-e\{k) 



^-^[^fe,4"|]- ^ 



■ CTfei, 



in which we have abbreviated the quantity 



(25a) 



(25b) 



The solution of this equation is found most easily using 
the time evolution operator e^^ht/h^ ^^^^ ^^ allow the 
electric field to have a small but finite frequency oj, taking 
the limit a; ^ at the end, 



Using Eqs. (11a) and (15) we can write down the equa- 
tion for Sf^n and equate the coefficients of crfe|| , reducing 

it to an equation for s^,, . It is important to note from 
Eq. (19a) that the coefficient of (Tfc|| contains only the 
Fourier components to = ±1, and therefore the sum in 
Eq. (15) will only contain Ci (besides Co)- This means 
that we have for s 
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Here the scattering time r is given by 
1 _ kn,W 

where the term in brackets can be expressed as 

^ 2 



Co - Ci = 




kp kp 



Using Eq. (9a) this gives us 5^11 as 
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As discussed above, unlike S^w given by Eq. (22), the 
perpendicular part Ski_ of the spin density matrix is in- 
dependent of the transition rate |C/fefe'p. 



C. Sfeii to zeroth order in scattering 



Finally, S'L is found from Eq. (lie) 

P\\JisS) = -P\\J(sl^l) 



(27) 



{21&J Since S'^j^ is known, we need to take the expression for 

Sf^j_, act on it with the scattering operator and project 
the resulting expression parallel to H. This will then be- 
come the source term for Sj^J. The details of this process 
(21b) arc given in Appendix B. The result is 



^(0) ^ eE-kXik) ^^ 



1 



2ft jj,w^o 2i{2vk — UJ — i-q) 

AO) 



<Jk\\- (28) 



This term is very similar to Sf,j_ and their averages over 
directions in momentum space are the same. 



D. Electrical conductivity 

Using the velocity operator in single-layer graphene, 
given by V ~ V (T, we can finally determine the electrical 
current. The current operator depends on the pseudospin 
and, following our reasoning so far, is decomposed into 
a parallel part (||) and a perpendicular part (_L). The 
expectation value of the current reads 



jx = -e hm 



_fk_ 
(27r)2 
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Sk\\ 
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(29) 



We convert the current tensor into the conductivity ten- 
sor <T using Ohm's law j = crE. The tensor a is diagonal, 
with a XX — <^yy and (Jxy = 0. The contribution to the 
conductivity due to S^ii (per valley and spin) is 



XX ^^ 



(30) 



^ord 

XX 



a^"^ behaves differently depending on the nature of scat- 



terers in the system.*^ For long range scatterers r = 
{4:h'^e'^vkF)/{niZ'^e*), with Z the atomic number and e 
the permittivity, thus CT°|jf oc (n/rii). Short range scat- 
terers give T = {2h'^e^vk^p)/{niZ'^e'^kF), but since in 
graphene hrr k kp as shown by Ando we still have 
the relationship a°'il^ oc {n/ui). 

The contribution due to Sj^_[ requires a careful evalu- 
ation of Ihnrj^i^^o J dk fo±/i'2vk — lo — irf). This is per- 
formed by replacing first 1 / {2vk — uj — irf) by m5{2vk — Lj) 
while keeping w 7^ 0, then evaluating the integral using 
the 5-function. (Otherwise one can obtain a negative 
conductivity.) We get 
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Note that the limits w — > and (3 ~ l/{k-QT) -^ 00 are 
not equivalent. Results of a similar or identical magni- 
tude have been found before in the absence of disorder 
and for n — only.'^®"^'* The contribution of S'L to the 

conductivity is equal to (per valley and spin) axx = cr^x ■ 
This result reflects the fact that Sj^l^ and S'L have the 
same angular average in momentum space. Remarkably 
this holds for a screened Coulomb potential regardless of 
the screening length. 



changing mixture of electrons and holes). Pseudospin 
non-conservation, which is crucial in determining the 
scattering-independent contribution to the conductiv- 
ity, occurs in graphene due to the mixing of electron 
and hole states contained in Hp, {electron-hole coher- 
ence). The steady state in graphene therefore involves 
non-conservation of the pseudospin due to the mixing 
of electron and hole states present in the Hamiltonian. 
The derivation shown in the preceding section demon- 
strates that the conserved pseudospin distribution gives 
rise to the fraction tT°|jf of the conductivity given by Eq. 
(30). cr°|jf corresponds to the ordinary electrical ("Boltz- 
mann") conductivity, which in the steady state is pro- 
portional to the carrier number density n and the char- 
acteristic scattering time r. The non-conserved pseu- 
dospin distribution gives rise to a second contribution 
to the electrical conductivity, cr'^x , which appears to 
be independent of the carrier number density and scat- 
tering time. However, our analysis shows that this is 
not the complete answer, since scattering between Ski_ 
and S'fcii produces an additional correction of order zero 

in the scattering potential, aJx which reinforces cr^x- 
Scattering from the non-conserved pseudospin distribu- 
tion into the conserved pseudospin distribution is repre- 
sented by P|[ J{Ski) in Eq. (lie). In that sense scatter- 
ing between the two distributions has a constructive ef- 
fect, and we emphasize that this reinforcement holds for 
screened Coulomb scattering, regardless of the screen- 
ing length (that is, both short-ranged and long-ranged 
impurity potentials.) Conversely, the additional correc- 
tion axx , which depends on the number density but not 
on the impurity density, is also a result of scattering. 
We emphasize that this reinforcement of the n- and r- 
independent contribution to the electrical conductivity in 
graphene in the weak momentum scattering limit was not 
found before and constitutes the main result of our work. 
This work demonstrates the unity behind two situations 
which until now appeared to be two different limits of 
the problem of determining the electrical conductivity of 
graphene. To date no approach has been put forward 
in which the Boltzmann and n- and r-independent con- 
tributions to the conductivity are treated on the same 
footing. 



VI. COMPARISON WITH SYSTEMS WITH 
SPIN-ORBIT INTERACTIONS 



V. DISCUSSION 

We would like to dwell on the role of the pseudospin 
in the electrical conductivity of graphene. In our anal- 
ysis the pseudospin density matrix Sk is quite gen- 
erally decomposed into a part Sk\\ representing con- 
served pseudospin (i.e., carriers which are either elec- 
trons or holes) and a part Sk± representing non- 
conserved pseudospin (carriers which are a continually 



The derivation of the steady-state density matrix in 
graphene presented in this work makes evident the many 
parallels which exist between the establishment of electri- 
cal currents in graphene and the establishment of steady- 
state spin densities and spin currents in semiconductors 
with strong spin-orbit interactions. In a recent paper^^ 
we discussed the nature of the steady state in systems 
with spin-orbit interactions, and showed that it is very 
different from the steady state established in usual charge 
systems. This difference is due to the presence of spin 



precession as a result of spin-orbit coupling. In the steady 
state in spin-orbit systems the spin density matrix is de- 
composed into a part representing conserved (i.e., not 
precessing) spin and a part representing precessing spin. 
The conserved spin distribution is responsible for the es- 
tablishment of a steady-state spin density, which is pro- 
portional to the carrier number density and the charac- 
teristic scattering time r. The precessing spin distribu- 
tion is responsible for steady state spin currents, which 
are independent of the scattering time and in two dimen- 
sions appear to be independent of the number density. 
Interestingly, the correction equivalent to SL vanishes 

in spin-orbit systems. Our recent work^® demonstrated 
that scattering from the conserved spin distribution into 
the precessing spin distribution produces a correction to 
the precessing spin distribution which in general acts to 
reduce spin currents, and in certain circumstances causes 
the spin current to be zero. In that sense scattering 
between the two distributions also has a destructive ef- 
fect. Furthermore, our work showed that this cancela- 
tion is due to the same physics that produces the vertex 
correction to the spin conductivity in Green's functions 
approaches. ^^ 

Interestingly, the reinforcement of cr^j. in graphcnc 
happens in a way that is very similar to the vanish- 
ing of the spin current mentioned above. In both cases 
the effect is caused by the term Py J{Sj,j_) appearing in 
Eq. (lie). Yet this term makes contributions of differ- 
ent signs in the two systems. We remark that Eq. (25) 
has been solved by applying the time evolution opera- 
tor. The product e~'^'=*/'^(7fc_Le*^'°*/'' is made up of two 
terms which have different dependencies on the angle of 
the wave vector k. The first term is proportional to (7k± 
and the second term is proportional to [Hk,(7k±]. The 
electrical current operator in graphcnc is isotropic so that 
only the first term contributes to the conductivity in this 
material. In contrast, the spin current operator is pro- 
portional to k. Thus only the second term contributes 
to the spin current in spin-orbit coupled systems. These 
results imply that the correction to the spin current is 
brought about scattering out of Sk\\ and into Skj_. On 
the other hand in graphene it is scattering out of Ski_ 
and into S^w that gives rise to the correction to electrical 
conductivity. From the discussion above we expect this 
reinforcement to appear as a result of the vertex correc- 
tion to the charge conductivity in graphene if a linear 
response approach based on Green's functions is used. 

In the field of spin transport a simple and elegant argu- 
ment has been formulated^'^ which explains why the spin 
current is necessarily zero in certain systems. Briefly, 
the spin current is proportional to the rate of change of 
one of the spin components, which must be zero in the 
steady state. A similar argument cannot be made for 
charge transport in graphene, where the charge current 
is proportional to the pseudospin and does not depend 
on the rate of change of any quantity in the steady state. 
This is a reassuring observation: if an analogous argu- 



ment could be made it would imply that the electrical 
conductivity of graphene vanishes identically, and this is 
evidently not the case. 

Finally, we would like to point out that in spin- 1/2 
electron systems with spin-orbit interactions, where the 
kinetic energy greatly exceeds the spin-orbit splitting, it 
is customary to treat the spin-orbit interaction as a per- 
turbation. Yet it is interesting to bear in mind that one 
expects results qualitatively similar to those in graphene 
for a Rashba-type Hamiltonian with a very large spin- 
orbit constant. 



VII. SUMMARY 

We have examined closely the nature of the steady 
state established in graphene in the presence of an elec- 
tric field. We have demonstrated that the steady state in 
this material contains important qualitative differences 
from the steady state in other known conductors. The 
principal reason behind this difference is the existence of 
a pseudospin degree of freedom, which is related to the 
coupling between electrons and holes contained in the 
Hamiltonian. In the weak momentum scattering regime 
there are two contributions to the electrical conductiv- 
ity in graphene, one linear in the carrier number den- 
sity and scattering time and one independent of both. 
These contributions can be identified with pseudospin 
conservation and non-conservation respectively, and are 
connected by scattering processes. Scattering between 
the non-conserved and the conserved pseudospin distri- 
butions doubles the contribution to the conductivity in- 
dependent of n and r. Moreover, the steady-state den- 
sity matrix in graphene displays remarkable similarities 
to the steady-state spin density matrix in systems with 
spin-orbit interactions. The contribution linear in n and 
T has an analog in the part of the spin density matrix 
which yields a steady state spin density, while the contri- 
bution independent of n and r is analogous to the part of 
the spin density matrix which yields a steady state spin 
current. The reinforcement of the n-and t- independent 
part of the conductivity is due to scattering between the 
conserved and non-conserved pseudospin distributions. 
This scattering also has an analogy in spin- 1/2 electron 
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ex- 



systems with spin-orbit interactions linear in k 
cept in those systems scattering between the conserved 
and non-conserved spin distributions causes the spin-Hall 
current to vanish. 
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We express pvi{r) in terms of its Fourier transform, using 
the convention (r|fc) = e"*'^'", as 

/Mw^/igj/H, «-•"-';'■' «*,n (A3) 

vi^here d — 2 for graphcnc. At this point we substitute the 
expression for the spatially-inhomogeneous electric field 
and focus on the part of the density matrix diagonal in 
vifave vector 



Tr(OS) = tr 



d'^k ^ ,, , eE dpa 



{2n) 
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h dk 
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This tells us that the source term in the kinetic equation 
is therefore {eE/K) ■ {dpo/dk). 



APPENDIX B: Sk\\ TO ZEROTH ORDER IN 
SCATTERING 

We will discuss in this section the details of the con- 
tribution S'j^ii , which is found from 
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APPENDIX A: STEADY STATE EXPECTATION 
VALUES 

Our aim is to derive a kinetic equation for the density 
matrix in the presence of an electric field E. In order 
to evaluate the electric field-induced source term in the 
kinetic equation we start from the Liouville equation with 
the perturbing Hamiltonian H^ = eE ■ r. The source 
term in the Liouville equation is 



ieE 



[r,Pa\ 



(Al) 



where po is the density operator in equilibrium, that is, 
in the absence of the external field. We are interested in 
the expectation value of an operator O (in our case the 
electrical current), which is found by taking the trace of 
this operator with the density matrix, and thus with the 
correction to the density matrix due to the electric field. 
This correction evidently depends on the source so that 
we consider first the trace of the operator O with the 
source term due to the electric field. This trace (denoted 
by Tr) is evaluated as follows 
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As mentioned in the main text, Sj^j_ is known. We act on 
it with the scattering operator and project the resulting 
expression parallel to H, yielding 
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The term ex vanishes in the angular integration, while 
the remaining term gives 



P\\J{Ski_) 



eE-kX{k) 



lim 



1 



(B3) 
This will then act as the source term for S'^J , and the 
equation is solved in the same way as in Sec. IV A 
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